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Abstract 

The dynamics of a class of cosmological models with collisionless mat- 
ter and four Killing vectors is studied in detail and compared with that of 
corresponding perfect fluid models. In many cases it is possible to identify 
asymptotic states of the spacetimes near the singularity or in a phase of 
unlimited expansion. Bianchi type II models show oscillatory behaviour 
near the initial singularity which is, however, simpler than that of the 
mixmaster model. 



1 Introduction 

In studies of the dynamics of spatially homogeneous cosmological models it is 
usual to choose a perfect fluid with linear equation of state to describe the 
matter. The book [[[9) provides an excellent guide to the subject. In view of 
the fact that this restriction is made so frequently in the literature, it is natural 
to pose the question to what extent the conclusions obtained would change if 
the matter model were chosen differently. In juj it was shown that in the case 
of collisionless matter described by the Vlasov equation significant changes can 
occur in comparison with the case of a perfect fluid. More specifically, it was 
shown that a solution of Bianchi type I exists whose qualitative behaviour near 
the initial singularity is different from that of any spacetime of that Bianchi type 
whose matter content is a fluid with a physically reasonable equation of state, 
linear or nonlinear. In the following this analysis will be generalized to show just 
how different models with collisionless matter can be from models with perfect 
fluid having the same symmetry. Differences are found in models of Bianchi 
type II (Theorem 4.2), Bianchi type III (Theorem 5.2) and Kantowski-Sachs 
models (Theorem 5.1). These concern both the initial singularity and phases 
of unlimited expansion. Perhaps the most striking case is that of the initial 
singularity in the Bianchi type II models, where we find persistent oscillatory 
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behaviour near the singularity. This is quite different from the known behaviour 
of type II perfect fluid models. 

Our results will also illuminate another matter. In [ pT| Lukash and Starobin- 
ski gave a heuristic analysis of a locally rotationally symmetric (LRS) model of 
Bianchi type I with collisionless matter consisting of massless particles. Their 
conclusion was that in the expanding direction the model would isotropize so 
that at large times it would look like a Friedman-Robertson- Walker model. On 
the one hand we are able to prove rigorously that the heuristic analysis of jll| 
gives the correct result. On the other hand we show that this result depends 
essentially on the assumption of a symmetry of Bianchi type I. If this symmetry 
type is replaced by Bianchi type II (keeping the LRS assumption and massless 
collisionless particles) then the anisotropy tends to a constant non-zero value at 
large times. 

The cosmological models studied in this paper are LRS spatially homoge- 
neous spacetimes with matter described by the Vlasov equation for massless 
particles. The reason for imposing the LRS condition is that it allows the Vlasov 
equation to be solved explicitly so that the Einstein- Vlasov equations reduce to 
a system of ordinary differential equations, albeit with coefficients which are not 
explicitly known and depend on the chosen initial data. The reason for choos- 
ing the particles to be massless is that this allows a reduction of the system of 
ODE similar to that carried out for perfect fluids with a linear equation of state 
by Wainwright and Hsu [ p0[ . It has not proved possible to analyse the global 
behaviour of solutions to our system of ODE completely. However a number of 
partial results have been obtained which show that there is considerable variety 
in the asymptotic behaviour of solutions near an initial singularity or during 
a phase of unlimited expansion. In particular, the reflection symmetric LRS 
Bianchi type I solutions with massless particles are analysed completely with 
respect to their asymptotic behaviour, thus improving markedly on the results 
obtained on that class of spacetimes in |l4| . 

The matter model used in the following will now be described. The matter 
consists of particles of zero rest mass which propagate through spacetime freely 
without collisions. Each particle is affected by the others only by the gravita- 
tional field which they generate collectively. The wordline of each particle is a 
null geodesic. Each geodesic has a natural lift to the tangent bundle of space- 
time. Thus the geodesic equation defines a flow on the tangent bundle. By 
means of the metric this may if desired be transported to the cotangent bundle 
and here it will be convenient to do so. The subset of the cotangent bundle 
consisting of all covectors obtained by lowering the index of future-pointing null 
vectors, which will be denoted by P, is invariant under the flow and thus the 
flow may be restricted to it. The basic matter field used to describe the colli- 
sionless particles is a non-negative real- valued function / on P which represents 
the density of particles with given position and momentum at a given time. 
Choosing appropriate coordinates x a on spacetime and letting (x a ,p a ) be the 
corresponding coordinates on the cotangent bundle, the manifold P can be co- 
ordinatized by (x a ,p a )- Here the convention is used that Greek and Roman 
indices run from to 3 and 1 to 3 respectively. We write t for x° and it is as- 



2 



sumed that t increases towards the future. The field equation for /, the Vlasov 
equation, says geometrically that / is constant along the geodesic flow. In the 
coordinates just introduced its explicit form is: 



where p° is to be determined from p a by the relation g a pp a p — and indices 
are raised and lowered using the spacetime metric g a p and its inverse. In order 
to couple the Vlasov equation to the Einstein equation, it is necessary to define 
the energy-momentum tensor. It is given by 



In fact for Bianchi models it is more useful to replace the coordinate components 
of the momentum used in these equations by components in a suitable frame. 
The only change in the equations is that the Christoffel symbols in the Vlasov 
equation are replaced by the connection coefficients in the given frame. For 
more information about the Vlasov equation in general relativity the reader is 
referred to || and 

Spatially homogeneous spacetimes fall into three broad classes, known as 
Bianchi class A, Bianchi class B and Kantowski-Sachs (see [[l9]). Each of the 
two Bianchi classes can be further divided into Bianchi types. A spatially ho- 
mogeneous spacetime in one of the Bianchi classes is called locally rotationally 
symmetric if it has, in addition to the three Killing vector fields needed for 
spatial homogeneity, a fourth one. This can only happen for certain symmetry 
types. In class A the Bianchi types which allow an LRS special case are I, II, 
VII , VIII and IX. In class B it is types III, V and VII^ which allow this @. The 
Kantowski-Sachs spacetimes automatically have a fourth Killing vector. There 
exist solutions of the Einstein- Vlasov equations with k = — 1 Robertson- Walker 
symmetry and these have, in particular, Bianchi type V and Bianchi type Vllh 
symmetry with any non-zero h. We did not attempt to ascertain whether there 
are other examples of solutions of these Bianchi types with LRS symmetry, and 
these types are not considered further in this paper. A spatially homogeneous 
solution of the Einstein- Vlasov equations has by definition the property that 
both the geometry and the phase space density of particles are invariant under 
the group action defining the symmetry type. A similar remark applies to an 
additional LRS symmetry. It would be nice if the invariance of / under the 
group in a Bianchi model could be expressed by the condition that / depends 
only on time and momentum when expressed with respect to a left-invariant 
frame on the group defining the symmetry. Unfortunately, as discussed in Q, 
this does not work in general. It does work for all LRS Bianchi models of class 
A and type III and for Kantowski-Sachs models |Q. This is the reason why 
LRS models are relatively tractable. In the following we consider LRS models 
which are of Kantowski-Sachs type, or of Bianchi type I, II, III, VIIq, VIII or 



In the next section it is shown how in the class of spacetimes of interest the 
Einstein- Vlasov equations with given initial data can be reduced to a system 



df/dt+ ( P a /p°)df/dx a + (r^p a py P )df/dp b = o 



(i) 





IX. 
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of ordinary differential equations. In fact two systems are needed. The first 
includes the solutions of types I, II, VII , VIII and IX while the second includes 
those of types I and III and the Kantowski-Sachs models. Note that the solutions 
of type I are represented in both systems and understanding the Bianchi I case 
is central to analysing the general case. The analysis of the Bianchi I system 
is carried out in the third section. This is then used in Sections 4 and 5 to 
obtain results on the first and second systems of ODE respectively. In the last 
section the results are summarized and their wider significance is examined. An 
appendix collects together some results from the theory of dynamical systems 
used in the body of the paper. 

2 Reduction to an ODE problem 

In a spacetime with Bianchi symmetry the metric can be written in the form 

ds 2 = -dt 2 +g ab {t)9 a ®9 h (3) 

where {0 a } is a left-invariant coframe on the Lie group G which defines the 
symmetry. The particular Bianchi type is determined by the structure constants 
of the Lie algebra of G. The extra symmetry which is present in the LRS case 
implies that the metric g a b(t) is diagonal, with two of the diagonal elements 
being equal [Q. Thus (||) simplifies to 

ds 2 = -dt 2 + a 2 (t)(^) 2 + b 2 (t)((9 2 ) 2 + (9 3 ) 2 ) (4) 

for two functions a(t) and b(t) of one variable. If k a is any Killing vector field 
then the function p a k a on the cotangent bundle is constant along geodesies and 
hence satisfies the Vlasov equation. Any function of quantities of this type for 
different Killing vectors also satisfies the Vlasov equation. The Killing vectors 
on a spacetime with Bianchi symmetry include those defined by right-invariant 
vector fields on the Lie group G but the result of evaluating a left-invariant one- 
form on one of these is not in general constant. Thus we cannot simply solve the 
Vlasov equation by choosing an arbitrary function of the components p a with 
respect to a left-invariant basis. However for the LRS spacetimes of Bianchi class 
A or type III considered here a function of the form f(t,pi,p2,P3) = fo{pi,p 2 + 
p 2 ) does satisfy the Vlasov equation and in fact is the most general solution 
with the full LRS symmetry [[l2f . Here p\ , P2 and P3 are the components of the 
momentum in the coframe {&"}■ Since / does not depend explicitly on time in 
this representation, the function /o can be identified with the initial datum for 
the solution of the Vlasov equation at a fixed time. A similar statement holds for 
Kantowski-Sachs spacetimes. The metric can be written in the form (|j) where 
9 1 is invariant under the symmetry group and 9 2 and 9 3 make up any (locally 
defined) orthonormal coframe on the two-sphere. The expression p\ + p 2 is not 
changed by a change in orthonormal coframe and so it makes sense to consider 
the above form of / in terms of fo in Kantowski-Sachs spacetimes as well. If / 
is of this form it satisfies the Vlasov equation. Thus the Vlasov equation has 
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been solved explicitly in the class of spacetimes to be studied. It remains to 
determine the form of the Einstein equations. In fact, one further restriction 
will be imposed. The distribution function given above is automatically an even 
function of p 2 and p% . However it need not be even in p\ . If it is even in p\ 



we say, as in 14 , that the solution is reflection symmetric. Only reflection 
symmetric solutions will be considered in the following. For convenience we say 
that a function of pi, p2 and p% which depends only on p\ and p\ +p| and which 
is even in p\ has special form. 

If the Einstein equations are split as usual into constraints and evolution 
equations then it turns out that in this class of spacetimes the momentum 
constraint is automatically satisfied. Only the Hamiltonian constraint and the 
evolution equations are left. The former is an algebraic relation between a, b 
and their time derivatives da/dt, db/dt. The latter provide ordinary differential 
equations for the evolution of a and 6 which are second order in time. It will 
be convenient to write these equations in terms of some alternative variables. 
Consider first the mean curvature of the homogeneous hypersurfaces: 

trjfe = -[a' 1 da/dt + 2b~ 1 db/dt] (5) 

A new time coordinate r can be defined by r(t) = — f t tik(t)dt for some arbi- 
trary fixed time to. In the following a dot over a quantity denotes its derivative 
with respect to r. Now define: 

q = b/a, 

m = -ex(a/6 2 )(trk)-\ 

N 2 = -e 2 a- 1 (trk)- 1 , (6) 

£+ = -3(6- 1 rf6/^)(trk)" 1 - 1, 

B = -b^itrk)- 1 

where e± and e 2 will be —1, or 1, depending on the symmetry type considered. 
The variables Ni, N2 and E + are closely related to the variables of the same 
names used by Wainwright and Hsu pCfl . (Note that we adopt the conventions 
of |2(J rather than those of [^9|, which differ by a factor of three in some places.) 

Two systems of ODE will now be considered, which between them are equiv- 
alent to the evolution part of the Einstein- Vlasov equations for all the relevant 
symmetry types. 

The first system is: 

q = Y, + q 

Nt = [_iJV 1 (JV 1 -4^ 3 ) + |(l-4S + + E5.)]JVi 

N2 = [-i7V 1 (7V 1 -47V 2 ) + i(l + 2E+ + S^)]7V 2 (7) 

±+ = |{|iV 1 2 + IiV 1 (iV 1 -4iV 2 )(l-2S + ) 

+ [-Mm - 4N 2 ) + i(l - S 2 + )][i(l - 2S+) - Q}} 
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Here Q is defined to be 



Q(q) = q 2 



pi) 1/2 dp 1 dp 2 dp 3 



ffo( Pi )(q 2 pl 



p\ 



-pl) 1/2 dp 1 dp 2 dp 3 



(8) 



where /o is a fixed smooth function of special form and compactly supported 
on R 3 . The Hamiltonian constraint is 



167rp/(tr£0 2 = -|iVi(JV x - AN 2 ) + |(1 - Ti%) 



(9) 



where p is the energy density and to take account of the positivity of p, only 
the region satisfying the inequality 

- |JV x (JVi - 4iV 2 ) + |(1 - £2_) > (10) 

is considered. Define submanifolds of this region by the following conditions: 



Si 


Ni 


= iV 2 = 





s 2 


JVi 


^0,7V 2 


= 


s 3 


Ni 


= 0,7V 2 




s 4 


iVk 




^0,N 2 


s 5 


iVi 


^0,N 2 


^0,N 2 



q 2 N x 



The submanifolds S\, S 2 , S3, S4 and S5 correspond to Bianchi types I, II, VIIo, 
VIII and IX respectively. To make the correspondence with spacetime quantities 
in these different cases (ei, e 2 ) should be chosen to be (0, 0), (1,0), (0, 1), (—1, 1) 
and (1, 1) respectively. Note that if q is replaced by q — q^ 1 in (0) an almost 
identical system is obtained, with the sign in the first equation being reversed. 
The second system is: 



q 

B 



\eB 2 



^(l-2£ + ) 2 ]B 



(11) 



|{-|eB 2 (l-2E + ) + [eS 2 



i(l-£ 2 )][i(l-2£ + )-Q]} 



where e belongs to the set {—1,0, 1}. Only the region satisfying the inequality 



2eB 2 



(1 - S 2 ,) > 



(12) 



is considered. The cases e = — 1, e = and e = 1 correspond to Bianchi type 
III, Bianchi type I and Kantowski-Sachs respectively. Note that the restriction 
of the system (uh to Si is identical to the system consisting of the first and third 
equations of (fnl) for e = 0. This restricted system will be referred to in the 
following as the Bianchi I system. It was introduced in section 6 of |l4| with 
slightly different variables. 

If a solution of (Q) and a fixed fo are given, it is possible to construct a 
spacetime as follows. Suppose that r = is contained in the domain of definition 
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of the solution. Since the system is autonomous this is no essential restriction. 
Choose a negative number Ho- Define 

p = (l/167r)flg[-|JVi(0)(JVi(0) - 4JV 2 (0)) + |(1 - E + (0))] (13) 

Let 

1= J fo(p l )[(q(0)) 2 pi+pi+pl} 1/2 dp 1 d P2 dp 3 (14) 

and define 

bo = p- x IW\q{0)) 1,A (15) 
In terms of these quantities we can define an initial metric by 

al(e x f + blue 2 ) 2 + (e 3 f) (is) 

Similarly, we can define an initial second fundamental form by 

- 1(1 - 2E + (0))if a^ 1 ) 2 + i ( i + X +{ Q))H bl((e 2 ? + (0 3 ) 2 ) (17) 

These data satisfy the constraints by construction. Consider now the spacetime 
which evolves from these initial data. It is of the form (||). For the Einstein- 
Vlasov system in a spacetime of the form (^|) with a fixed time-independent 
distribution function is a system of second order ODE which has solutions cor- 
responding to data for (a,b,da/dt,db/dt). These data can be chosen so as 
to reproduce the data of interest for the Einstein- Vlasov system by choosing 
da/dt = i(l - 2S+(0)) J ff o a and db/dt = -|(1 + Y, + (0))H o b for t = t . This 
spacetime defines a solution of (0) via (||). (Note that t = to corresponds to 
t = 0.) Thus the two solutions are identical. In this way a spacetime has been 
constructed which gives rise to the solution of (Q) we started with. This space- 
time may be obtained more explicitly if desired. In order to do this, first solve 
the equation: 

«9 T (trfc) = _[-|JVi(JVi - 4N 2 ) + |(2 + E + )]trfc (18) 

with initial data Ho- Then p can be obtained from the Hamiltonian constraint 
. The definition of p in terms of /o can then be combined with q to give a 
and b as in (p5|). Finally t can be obtained from tr/c. All the considerations here 
in the case of (Q) are equally applicable in the case of (|i~l|). The analogue of 
equation (|TJ) is 

dr(trk) = -[eB 2 + |(2 + E + )]tr* (19) 

Solutions of the Einstein equations with matter described by a perfect fluid 
with a linear equation of state p — (j—l)p which belong to one of the symmetry 
types studied in the case of collisionlcss matter in the following can be described 
by equations very similar to (0) and (pr|). The similarity is particularly great 
in the case 7 = | (radiation fluid). In that case the only difference is that 
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the function Q{q) should be replaced by the constant value |. This leads to a 
decoupling of the first equation in each system, so that it is possible to restrict 
attention to the remaining equations when investigating the dynamics. (This 
last remark also applies to the system obtained for other values of 7.) The 
equation for q can be integrated afterwards if desired. 

In [H it was proved that Q{q) as defined in (j^) tends to zero as q tends to 
zero and that if Q(0) is defined to be zero the resulting extension of Q is C 1 with 
Q'(0) = 0. This means in particular that the dynamical system (0) has a well- 
defined C 1 extension to q = 0. In a similar way it can be shown that if a function 
Q is defined by Q(q) = Q(q) then Q can be extended in a C 1 manner to q = 
in such a way that Q(0) = 1 and Q'(0) = 0. For 1 - Q = (p - TjVp = 2T$/p 
and this last expression is 0(q 4 ^ 3 ) as q — ► by Lemma 4.2 of [H. By using 
a coordinate q — q/(q + 1) it is possible to map the system (fjj) with q ranging 
from zero to infinity onto a region with q ranging from zero to one. Moreover 
the system extends in a C 1 manner to the boundary components q = and 
q = 1. The coordinate q has been introduced purely to demonstrate that the 
system ([?]) can be smoothly compactified in the g-direction. For computations 
it is more practical to use the local coordinates q and q. In particular, these 
considerations allow us to regard the Bianchi I system as being defined on a 
compact set. 

The compactification of the system (Q) is defined on a region with boundary. 
Different parts of the boundary are given by q = 0, q = 1 and the case of equality 
in (10). The complement of the boundary will be called the interior region in the 
following. A solution which lies in the interior region corresponds to a smooth 
non- vacuum solution of the Einstein- Vlasov equations. A solution which lies 
in the part of the boundary where ( |l0| ) becomes an equality corresponds to a 
solution of the vacuum Einstein equations. A solution which lies in the part of 
the boundary given by q = or q = 1 corresponds to a distributional solution of 
the Einstein- Vlasov equations, as will be explained in more detail below. The 
system (^Tj) can be compactified in a way very similar to that taken in the case 
of (0). The comments on the interpretation of different types of solutions of 
the compactification of (0) just made also apply to the compactification of (|TT|), 
with ( |lG| ) being replaced by jl^). 

Consider now the stationary points of the system (Q), or rather of its com- 
pactification. (This distinction will not always be made explicitly in what fol- 
lows.) In section 4 it will be shown that all stationary points where q has a 
finite non-zero value belong to the subset S\ corresponding to solutions of type 
I. In particular, they correspond to stationary points of the Bianchi I system, 
which will be studied in detail in the next section. 



3 The Bianchi I system 

It turns out that the Bianchi I system plays a central role in the dynamics of 
solutions of the systems (0) and (|ll|). In this section the asymptotic behaviour 
of solutions of this system is determined, both for r — ► — 00 (approach to the 
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Figure 1: The (g, E + ) plane and the fixed points for Bianchi type I. 

singularity) and for r — > oo (unlimited expansion). 

The first step in analysing the Bianchi I system is to determine the stationary 
points. This can be done using the fact, proved in O], that Q is strictly 
monotone for q > so that there is a unique qo with Q(qo) = ■§■ With this 
information it is straightforward to show that the coordinates of the stationary 
points in the (q, £+) plane are (go, 0), (0,-1), (0, ^), (0, 1), (oo, —1) and (oo, 1). 
Here q = oo is to be interpreted as q = or q = 1. Call these points P\,.. . ,P§ 
respectively (see figure 1). The next step is to linearize the system about the 
stationary points. Recall that a stationary point is called hyperbolic if none of 
its eigenvalues are purely imaginary. In the following we call a stationary point 
degenerate if it is not hyperbolic. The point P\ is a hyperbolic sink while P4 is 
a hyperbolic source. The points P2, P3 and Pq are hyperbolic saddles while P5 
is degenerate, with one zero eigenvalue. 

Before proceeding further, we state the main result of this section. 

Theorem 3.1 If a smooth non- vacuum reflection-symmetric LRS solution of 
Bianchi type I of the Einstein- Vlasov equations for massless particles is repre- 
sented as a solution of ([?]) with N\ = N2 = then for r — > 00 it converges to 
the point P\. For r — > — 00 either 

(i) it converges to P\ and in that case it stays for all time at the point Pi or 

(ii) it converges to the point P 3 and it belongs to the unstable manifold of P 3 
or 
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(iii) it converges to P4 

All of these cases occur, and (iii) is the generic case in the sense that it occurs 
for an open dense set of initial data. 

This will be proved in a series of lemmas. Terminology from the theory of 
dynamical systems which may be unfamiliar to the reader is explained in the 
appendix. 

Lemma 3.1 If a solution of the Bianchi I system in the interior enters the region 
£+ > 1/2 then for r — > —00 it belongs to case (iii) of Theorem 3.1. A solution 
of the Bianchi I system in the interior has no w-limit points with S + > 1/2. 
Proof A solution of the Bianchi I system satisfies £ < — — £+)<? when S + > 
1/2 and so for any solution which enters the given region, S + is nondccrcasing 
towards the past and it is in the region for all earlier times. If S + did not 
tend to 1 as r — > —00 then we would have S + < — C < at early times, a 
contradiction. Once we know that X+ — > 1 as r — > —00 it follows immediately 
that q — * 0. Thus the solution converges to P4. Consider now the forward time 
direction. Since is positive, q is increasing. This means that Q is increasing. 
The inequality X < — CQ for a constant C > then shows that the solution 
must leave the region of interest in finite time, so that there can be no w-limit 
point with S + > 1/2. 

For convenience an interior solution which does not tend to P4 as r — > — 00 will 
be called exceptional. Thus Lemma 3.1 says that an exceptional solution cannot 
intersect the region S + > 1/2. 

Lemma 3.2 The a- limit set of an exceptional solution cannot intersect the 
boundary at any point except P3. If it does intersect the boundary at P 3 it 
belongs to case (ii) of Theorem 3.1 as r — > —00. The w-limit set of any interior 
solution cannot intersect the boundary at all. 

Proof Let (q, S + ) be a point of the a- limit set of an exceptional solution which 
lies on the boundary. If q = 00 then the whole orbit passing through that point 
belongs to the a-limit set. This implies that the solution must intersect the 
region E + > 1/2, a contradiction. Thus in fact q < 00. If S + = —1 then all 
points with S + = — 1 must be in the a- limit set, in particular P 5 . By Lemma 
A2 of the appendix, it follows that a point of the centre manifold of P5 lies in 
the a-limit set. However, this centre manifold is given by q = 00 and so we 
again obtain a contradiction. Hence X+ > — 1. If q = and S + < 1/2 then all 
points satisfying these conditions must be in the a-limit set, in particular P2. 
But then an application of Lemma Al of the appendix leads to a contradiction. 
Thus no point on the boundary other than P 3 is possible. A further application 
of Lemma Al shows that in this case the solution must lie on the unstable 
manifold of P 3 . In a similar way it is possible to show that if any point of the 
boundary belonged to the w-limit set of an interior solution then some point 
with S + > 1/2 would do so. However we know from Lemma 3.1 that this is 
impossible. 

Proof of Theorem 3.1 First the Poincarc-Bendixson theorem will be applied 
to the restriction of the Bianchi I system to the interior with the point Pi 
removed. In general the a- and w-limit sets of an orbit of a dynamical system 
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can be very complicated, but in two dimensions (and the Bianchi I system is two- 
dimensional) things are a lot simpler. Complicated situations are still possible 
and these play an important role in Hilbert's sixteenth problem (see e.g. 
p. 104). However many pathologies are ruled out by the Poincare-Bendixson 
theorem, which is stated in the Appendix (Theorem A2). 

Given an interior solution, suppose that P\ does not belong to the w-limit 
set. By Lemma 3.2 no point of the boundary belongs to the w-limit set either. 
Since Pi is a hyperbolic sink it follows that there must be a neighbourhood of Pi 
which does not intersect the w-limit set. Thus the solution remains in a compact 
set of the interior with the point P\ removed as r — > oo. Then Theorem A2 
implies the existence of a non-stationary periodic orbit of the Bianchi I system. 
In fact the existence of periodic solutions of the Bianchi I system can be ruled 
out by the presence of a Dulac function. (For a discussion of this concept see 
p9|.) Define a function F(q) by 



and denote the vector field defining the Bianchi I system by X. Then div(GA) 
is negative. In fact it is a constant multiple of g _1 F 1 / 2 (l — £j_) _1 / 2 (2 — £+). 
This means that G is a Dulac function. It follows that the Bianchi I system 
has no periodic solutions. It can be concluded that Pi does lie in the w-limit 
set. But since Pi is a hyperbolic sink, this implies, via the Hartman-Grobman 
theorem (cf. Theorem Al), that the w-limit set consists of Pi alone, which 
proves the first part of the theorem. 

To prove the remainder of the theorem we can assume without loss of gen- 
erality that the solution is exceptional and that it does not lie on the unstable 
manifold of P3. If Pi were not in the a-limit set then we would get a contra- 
diction by the Poincare-Bendixson theorem and the absence of periodic orbits. 
Hence Pi must belong to the a-limit set, and since Pi is a hyperbolic sink, the 
only possibility left is case (i) of the theorem. 

The conclusion of this theorem can be summarized in words as follows. All solu- 
tions isotropize in the expanding direction. The initial singularity is gcncrically 
a cigar singularity but there are exceptional cases where it is a barrel or point 
singularity. (For this terminology see p9[ , p. 30.) Note for comparison that 
if the Vlasov equation is replaced by the Euler equation for a fluid satisfying 
a physically reasonable equation of state then there are no barrel singularities 
and all solutions which are not isotropic have cigar or pancake singularities (see 
[pT[). The pancake singularities are as common as the cigar singularities. In 
particular this means that for fluid solutions of Bianchi type I cigar singulari- 
ties are not generic. All solutions isotropize in the expanding direction. Note 
that the 'reasonable' equations of state include those of the form p = kp with 




(20) 



Then Q = (q/F)F'. For q > and |S+| < 1 let 



G(g,S + )=q- 1 P 1 / 2 (l-S 2 h )- 3 / 2 



(21) 
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< k < 1. The solutions of the Einstein- Vlasov equations approach an isotropic 
fluid solution with equation of state p = ^p in the sense that the tracefree part 
of the spatial projection Ty of the energy-momentum tensor divided by the en- 
ergy density p approaches zero, while trT/p = 1. The latter relation is always 
true for kinetic theory with massless particles and for a radiation fluid (equation 
of state p = 

4 Other class A models 

This section is concerned with the models of class A, as described by the system 

(0) . Only limited statements will be made about types VIII and IX. Even in the 
a priori simpler case of a perfect fluid with linear equation of state it is difficult 
to analyse LRS models of type VIII and IX. (For information on what is known 
about that case, see 

@, EI and section 8.5 of |§). A 

major difficulty is that 

in these cases the domain of definition of the dynamical system is non-compact. 
This allows the possibility that there may be anomalous solutions similar to 
those encountered in . Type I was analysed in the previous section and we 
will see that the analysis of type VIIo can be reduced to that case in a relatively 
straightforward way. The most interesting results are obtained for Bianchi type 
II. 

We start with a theorem on Bianchi type VIIo, which is a close analogue of 
Theorem 3.1. 

Theorem 4.1 If a smooth non- vacuum reflection-symmetric LRS solution of 
Bianchi type VIIo of the Einstein- Vlasov equations for massless particles is 
represented as a solution of (fj]) with iVi = then for r — > oo the pair (q, E + ) 
converges to (qo,0) while N 2 increases without limit. N 2 tends to zero as r — * 
—00 while the pair (q, S + ) cither 

(1) converges to Pi and in that case it stays for all time at the point Pi or 

(ii) converges to the point P3 and belongs to the unstable manifold of P3 or 

(iii) converges to P4 

All of these cases occur, and (iii) is the generic case in the sense that it occurs 

for an open dense set of initial data. 

Proof When Ni = the third equation in ([?]) becomes 

N 2 = i(l + S+) 2 A^ 2 (22) 

while the equations for E + and q do not involve N 2 . The latter equations form 
a subsystem which is identical to the equations for Bianchi type I, so that the 
situation is again as in figure 1. The qualitative behaviour of their solutions has 
been analysed in Theorem 3.1. All that remains to be done is then to put that 
information into equation (^2|) and read off the behaviour of N 2 . The expression 
(1 + E+) 2 is strictly positive for a non- vacuum solution of type VIIo, due to (|l0|). 
Moreover it is bounded by four. Thus the solution has the property that the 
sign of N 2 remains constant and the solution exists globally in r. It is also clear 
that N 2 — ► 00 as t — ► 00 and that N 2 — > as r — > —00. 



12 




Figure 2: The (q, E+,i\Ti) space and the fixed points for Bianchi type II; the 
lower two figures show the phase portraits on the end-faces q = 0, 1 

There is a simple explanation for the close relation between the Bianchi I and 
Bianchi VIIo solutions. They arc in fact the same spacetimes parametrized in 
two different ways. The full four-dimensional isomctry group has a subgroup of 
Bianchi type I and a one-parameter family of subgroups of Bianchi type VIIo. 

Next we turn to the solutions of type II. It will be shown that the stationary 
points of (the compactification of) (0) which lie in the closure of S2 are the points 
Pi, . . . , Pq which we know already together with one additional point P7, which 
has coordinates (00, |) (see figure 2). The corresponding distributional 
solution of the Einstein- Vlasov equations will be discussed in detail below. 

Theorem 4.2 If a smooth non-vacuum reflection symmetric LRS solution of 
Bianchi type II of the Einstein- Vlasov equations for massless particles is repre- 
sented as a solution of (0) with 7V2 = then for t — > 00 the solution converges 
to P7. For t — ► — 00 either: 

(i) the solution converges to Pi or 

(ii) the a-limit set of the solution consists of the points P2, Pi, P5 and P§ 
together with the orbits connecting Pi to P5, P5 to Pq and Pq to P4 in the 
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set iVi = and the stable manifold of P4 which connects P4 with P2 via the 
vacuum boundary. In particular nminf T — >--oo S + = — 1, limsup,.^.^ £+ = +1, 
liminfT-^-oo q(r) — —00 and limsup T _ > _ 00 q(r) — 00. 
Both of these cases occur and (ii) is the generic case. 

This theorem shows that while models of Bianchi type II have simple behaviour 
in the expanding phase, all tending to a single attractor, the behaviour near the 
initial singularity is in general oscillatory and quite different from the Bianchi 
type I case. Note also that the models of type II do not isotropize as r — * 00, 
which is another important difference from the type I models. 

A first important step in proving Theorem 4.2 is to use the identity 

d/dr^Nx) = q^Nxi-lN^Nx - 4A 2 ) + ±(1 - E* ) + fE^] (23) 

which holds for any solution of (0). This is done in the following lemma. 
Lemma 4.1 For any solution in the interior of S2 the following statements hold. 
As t — > 00 q tends to 00. Either the a-limit set is contained in the set q = or 
it contains one of the points Pi , . . . , P 6 ■ 

Proof For a non-vacuum solution strict inequality holds in ([l(]) and hence 
q 4 / 3 Ni is strictly increasing where it is non-zero. This means that as long as q 
is finite and the Bianchi type is II this quantity is always increasing. Since £2 
is compact, solutions exist globally in r. As r tends to plus or minus infinity 
the solution must go to the boundary of S^. Equation ( p3| ) shows that if q i ^N 1 
tends to a finite non-zero limit in either time direction then El is integrable 
on a half-infinite time interval. The derivative of this quantity is bounded and 
these two facts together imply that it must tend to zero in the limit. The same 



argument applies to the quantity appearing in (10) and so it must also tend to 
zero in the limit. Under these conditions N\ — > and E+ — » |r, a contradiction. 

It can be concluded that lim T ^ 0O (g 4 / 3 A''i)(r) = 00 and lim T ^_ 00 (Q ,4 / 3 A^i)(r) = 
0. From the first of these statements and the boundedness of N\ it follows 
that q tends to 00 as t tends to 00. The fact that q 4 ^ 3 Ni tends to zero in the 
contracting direction implies that the a-limit set is contained in the union of the 
sets q — and Ni = 0. Suppose the a-limit set contains some point for which 
q 7^ 0. This must belong to the Bianchi I set. Thus the a-limit set contains a 
solution of Bianchi type I. Using Theorem 3.1, we conclude that the a-limit set 
contains one of the points P±, . . . , Pq. 

The next lemma gives information about the nature of the stationary points 
on 52- We already know from Lemma 4.1 that these stationary points can only 
occur for Ni = 0, q = or q = 00. The stationary points Pi, . . . , Pq will be 
investigated first. The equations for q = and q — 00 will be studied in detail 
later. 

Lemma 4.2 The stationary points Pi, . . . , P4 and Pq of the restriction of the 
system (0) to S2 are hyperbolic saddles, while P 5 is degenerate. The stable 
manifold of Pi is given by = 0. The stable and unstable manifolds of P 2 are 
given by S + = —1, Ai =0 and q = respectively. The stable manifold of P3 
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is given by q = 0. The stable and unstable manifolds of P4 are given by q = 0, 
Nf = |(1 — E^) and N% = respectively. The stable and unstable manifolds 
of P 6 are given by iV 2 = ~(1 — E^) and and 9 = 00, Ni — respectively. The 
unstable manifold of P5 is given by iV 2 = 4(1 — £+). The set g = 00, ATi = is 
a centre manifold for P 5 . 

Proof All that needs to be done is to compute the linearizations of the system 
about the given points and to note that the manifolds named in the statement 
of the theorem are all invariant. Linearizing the restriction of (fj]) to JV2 = 0, 
and setting Ni = in the result, gives the system (a bar denotes a linearized 
quantity) : 

dq/dr = Y. + q + qT, + 
dNx/dr = i(l-4£ + + £2_)jVi (24) 

dt + /dr = [-±(l + £ + -3£ 2 + ) + £ + Q(< Z )]£+-±Q'( (? )(l-£ 2 h )g 

The linearization about P\ has eigenvalues -| and — | ± \ §<2'oQ'(So)- The 
invariant subspace of the linearization corresponding to the eigenvalues with 
negative real parts is the tangent space to N\ = 0. The linearizations about 
P-2, P 3 and P 4 are diagonal with diagonal entries (—1,2,1), {\,— \, — \) and 
(1,— |, i) respectively. Since P5 and P 6 lie at q — 00, we must change to the 
coordinate q to study the linearizations at these points. They are diagonal with 
diagonal elements (1,2,0) and (— 1,— §,§)■ 

Next the limiting systems for q — and q — 00 will be examined. 
Lemma 4.3 Consider the restriction of the system (^) to the set given by the 
equations ^2=^ = 0. For any solution which does not belong to the vacuum 
boundary and which does not satisfy N\ = 0, the a-limit set is the point P2 and 
the w-limit set is the point P3. 

Proof First it will be shown that the solution cannot be stationary. The equa- 
tion for E + shows that S+ > if S+ < |. Thus at a stationary point E + > |. 
On the other hand, the equation for N± shows that at a stationary point 

(£+ - 2) 2 = 3 + Nl > 3 (25) 

Using the fact that S + < 1 it follows that S + < 2 — V3 < i. Next, it follows 
from Theorem 3.1 on p. 150 of Q| that the solution cannot be periodic. It 
can be concluded using the Poincare-Bendixson theorem (Theorem A2 of the 
appendix) that the a- and cj-limit sets are contained in the boundary of the 
region. The behaviour of solutions on the boundary is easily determined. The 
nature of the stationary points on the boundary can be read off from Lemma 
4.2. P 2 is a hyperbolic source, P 3 is a hyperbolic sink and P4 is a hyperbolic 
saddle. The last fact means, using Lemma Al of the appendix, that P4 cannot 
be in the a- or w-limit set unless P2 or P3 is also. Thus it can be concluded that 
the a- and w-limit sets must contain either P2 or P3 and then the conclusion 
follows easily. 

The system given by N2 — and q — 00 is more complicated. 
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Lemma 4.4 Consider the restriction of the system (m) to the set given by the 
equations -/V2 = and q = 00. For any solution which does not belong to the 
vacuum boundary and which does not satisfy Ni — 0, the a-limit set consists 
of all points which are either on the vacuum boundary or satisfy Ni = 0. The 
w-limit set is the point P7. 
Proof Define a function Z by 

Z = Nl' 2 [-±Nt + |(1 - S 2 + )]3/4(i _ (26) 

This function is well-defined and continuous on S2 and smooth away from Ni = 
and -5-/V1 + |(1 - £+) = 0. Its derivative is given by 

d T Z = - iS + )- 1 [|(5S + - l) 2 + (-37V 2 + E 2 . - 15E+ + 4)(1 - Q)] (27) 

The restriction of Z to the set q = 00 is non-decreasing along solutions as a 
consequence of (|27|). Moreover, it is strictly increasing unless £+ = g. When 

£+ = i it follows from (0) that £ + ^ unless iVi = ^p. Thus apart from the 

stationary solution at the point P7 with coordinates (00, 2 -^p, |), the function 
Z is strictly increasing along any solution with q = 00. It follows that P7 is the 
w-limit point of all solutions. The function Z attains its minimum precisely on 
the boundary of the region where the system is defined and hence the a-limit set 
of any solution is contained in this boundary. The only stationary points on the 
boundary are P5 and P&. From Lemma 4.2 it follows that both are saddle points 
of this system. (P5 is degenerate while Pq is non-degenerate.) This suffices to 
show, using Lemma Al and Lemma A2 of the appendix, that the a-limit set 
consists of the entire boundary. 

Proof of Theorem 4.2 By Lemma 4.1 the oj-limit set of any solution consists 
of points with q — 00. Then Lemma 4.4 shows that either P7 belongs to the 
w-limit set or that the w-limit set consists entirely of points with q — 00 and 
TVi = or — iiV 2 + !(l — £*]_) = 0. A calculation of the linearization of (Q) around 
P7 shows that this point is a hyperbolic sink. Hence if P7 belongs to the w-limit 
set this set consists of P 7 alone. It remains to rule out the other possibility 
where the solution has an w-limit point on the boundary of the intersection of 
S2 with q = 00. In that case Lemma Al applied to the point Pq and Lemma 
A2 applied to the point P5 show that the w-limit set contains the whole of this 
boundary. It will now be shown using (^) that this leads to a contradiction. 
There exist <5i > and M > such that if |E + — g| > 61 and q > M the 
right hand side of ( p7[ ) is positive. This is because the first term dominates the 
second. By reducing Si and increasing M if necessary it can be ensured that 
there exist positive constants 771, 1J2 and 82 such that Z~ x d T Z can be bounded 
below by 771 as long as |£+ — i| > Si and q > M and E+ > 772 for |£+ — -|| < 6±, 

\Ni - > 8 2 and q > M. Finally, given 773 > there exists S 3 > so that 
E < 773 for |E| > 1 — ^3 and q > M. At sufficiently late times the solution lies in 
the region q > M. Moreover, under the present assumption on the w-limit set it 
cannot enter the neighbourhood of P7 defined by Si and 82 ■ Each time it crosses 
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the strip defined by |S+ — g| < Si at a sufficiently late time it must enter the 
region E + > I — S3 before it can return to the strip. It must spend a long time in 
the region E + > 1 — S3 (due to the smallness of 773). This time can be bounded 
below by Cr]^ 1 for a constant C > 0. During that time \ogZ must increase 
by at least (7(771/773). On the other hand \ogZ can only decrease while it is in 
the strip. It stays there for a time at most 5i/r]2 and can decrease by at most 
CSxjrfz- Thus the net change of Z for each time it enters the strip is at least 
C( r ?i/ r ?3 — S\/r}2)- If 773 is chosen small enough this will be bounded below by a 
positive quantity. Since the solution must, under the given assumptions, enter 
the strip infinitely often, this gives a contradiction. The proof of the statement 
about the w-limit set is now complete. 

Suppose that the a-limit set contains a point with q = and N\ 7^ 0. Then 
by Lemma 4.3 it contains P2 and P3 or P4. Lemma 4.1 then shows that at least 
one of Pi, . . . ,Pq is contained in the a-limit set. If the a-limit set contains Pi 
then cither the solution lies in the unstable manifold of Pi, which gives case 
(i) of the theorem, or the a-limit set contains points on that unstable manifold 
other than Pi itself. But since these satisfy neither Ni = or q = this is 
a contradiction. If the a-limit set contained points with Ni = with q finite 
and < 1 it would contain Pi, leading once again to a contradiction. If it 
contains P4 it follows from Lemma Al and what has just been said that it must 
contain either P3 or Pq. It must also contain P2. However, if it contained P3 it 
would, by another application of the same lemma contain Pi, a contradiction. 
On the other hand, if it contains Pi it must contain P4 and P5. If it contains 
P5 it must contain Pq and vice versa, by Lemma Al and Lemma A2. On the 
other hand, Lemma Al shows that if the a-limit set contains P5 or Pq it must 
contain P2 or P4. It also follows from these applications of the lemmas of the 
appendix that the relevant connecting orbits are contained in the a-limit sets. 

Now a spacetime corresponding to the point P7 will be determined (in figure 

2, this space-time follows a straight line at constant (£+,iVi) into P7.). From 

3 

equation ( |l9| ) it follows that trfc = H a e~ 5 T . Putting this in the equation relating 
t and r shows that trk = — §i _1 . Putting this in the third equation of (|^) gives 
b = b^t 2 / 3 . Equation (|^) implies that a = dot 1 / 3 . Finally, the second equation 
of (|6|) leads to the relation a = ^0-b 2 . Choosing an explicit representation of 
a Bianchi type II frame leads to the metric: 

ds 2 = ~dt 2 + \B 2 t 2 l 3 (dx + zdy) 2 + Bt 4 / 3 {dy 2 + dz 2 ) (28) 

where B is a constant. This metric is invariant under the homothety t 1— > At, 
x 1 ► A 2 / 3 x, y i-» A 1 / 3 y, z i~> A 1 / 3 z. It follows that td/dt + ^xd/dx+ \{yd/dy + 
zd/dz) is a homothetic vector field and that this metric is self-similar. It 
satisfies the Einstein equations with an energy-momentum tensor whose only 
non- vanishing components are p and T\\. These two are equal and are pro- 
portional to t~ 2 . This can be interpreted as a distributional solution of the 
Einstein- Vlasov equations with massless particles where the distribution func- 
tion is of the form f(pi,P2,P3) = fi{pi)S{p2)S{p3,) . (Note that a distributional 
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/ of this kind defines a dynamical system just as a smooth / does so that the 
solution can be represented in figure 2.) The exact form of the function /i is 
unimportant. Only the integrals J fi(j?i)pidpi and J fi{pi)Pidpi influence the 
energy-momentum tensor. Related to this fact is that the same spacetime can 
be interpreted as a solution of the Einstein equations coupled to two streams 
of null dust moving in opposite senses in the x 1 -direction. This corresponds to 
choosing f\ = (S(pi) + 8(—pi)) instead of a smooth function. The sum of two 
Dirac measures is necessary to preserve the reflection symmetry. This space- 
time has previously been considered by Dunn and Tupper|| in the context of 
cosmological models with electromagnetic fields, although it had to be rejected 
for their purposes since no consistent electromagnetic field existed. 

The monotone function Z which plays a crucial role in the proof of Theorem 
4.2 is rather complicated and so is unlikely to be found by trial and error. We 
found it by means of a Hamiltonian formulation of the equations for q = oo. 
Once the function was found for q = oo it was extended so as to be indepen- 
dent of q. In developing the Hamiltonian formulation we followed the treatment 
of Uggla in chapter 10 of A key point is that the energy density of a 

distributional solution of the Einstein- Vlasov system where the pressure is con- 
centrated in one direction can be related in a simple way to q. The function Z is 
the Hamiltonian for the (time dependent) Hamiltonian system. It was also the 
construction of Z which led us to discover the self-similar solution corresponding 
to the point P7. 

The picture obtained in Theorem 4.2 is quite different from that seen in LRS 
Bianchi type II solutions with a perfect fluid with linear equation of state as 
matter model (see [ pj| , chapter 6). There generic solutions are approximated 
near the singularity by a vacuum solution (the type II NUT solution) and there 
is no oscillatory behaviour. In the expanding direction the fluid solutions are 
also all asymptotic to a self-similar solution (the Collins-Stewart solution) but 
this solution has a different ratio of shear to expansion than the solution cor- 
responding to the point P 7 . Moreover the pressure is highly anisotropic in the 
latter solution. 

5 Kantowski-Sachs and Bianchi type III models 

In this section information will be obtained on Kantowski-Sachs models and 
models of Bianchi type III which is as complete as that obtained on models of 
Bianchi type I in section 3. 

Theorem 5.1 If a smooth non- vacuum reflection symmetric Kantowski-Sachs 
type solution of the Einstein- Vlasov equations for massless particles is repre- 
sented as a solution of ( |Ti"| ) with e = 1 then for r — > —00 either 

(i) it converges to Pi 

(ii) it converges to the point P3 and it belongs to the unstable manifold of P3 
or 

(iii) it converges to P4 
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All of these cases occur, and (iii) is the generic case in the sense that it occurs 
for an open dense set of initial data. 

Proof The inequality d T B > jB shows that B decreases towards the past. It 
follows that as r decreases the solution remains in a compact set and hence 
that the solution exists for all sufficiently negative r. Using the inequality again 
shows that B — > exponentially as r — > -co and the a-limit set lies in the set 
B = 0. The latter can be identified with the Bianchi I system. The a-limit set 
contains the image of a solution of the Bianchi I system and hence, by Theorem 
3.1 contains either Pi or some point of the boundary of the Bianchi I system. 
Each of the stationary points Pi , . . . , P§ , considered as stationary points of 
( pr| ) , has a linearization which differs from its linearization within the Bianchi I 
system by the addition of an extra eigenvector with a positive eigenvalue. It can 
be concluded from this that P4 is a hyperbolic source. Moreover, by Lemma Al 
and Lemma A2, if any point of the boundary other than P3 lies in the a-limit 
set, then P4 must also lie in the a-limit set. Hence in this case the a-limit set 
consists of P4 alone. Moreover, if P3 lies in the a-limit set then the solution must 
lie on its unstable manifold. The only remaining possibility is that the a-limit 
set consists of Pi alone, and that the solution lies on the unstable manifold of 
Pi- 

No statement is made here about the behaviour as r — > 00. In fact any solution 
of (^) with e = 1 tends to infinity in finite time. However this is not a problem 
from the point of view of understanding the spacetime. It is known that the 
Kantowski-Sachs models recollapse [|J. Thus there is no infinitely expanding 
phase and a final singularity looks like the time reverse of an initial singularity. 
One interesting question which we do not attempt to tackle here is whether 
there is an interesting correlation between the behaviour near the initial and final 
singularities. For each individual singularity the picture is essentially identical to 
that seen in the singularity of Bianchi I models. The system for a radiation fluid 
can be analysed in the same way, reducing the dynamics near the singularity to 
that of the corresponding Bianchi I system. The differences between radiation 
fluid and kinetic models are similar in both cases. 

Theorem 5.2 If a smooth non- vacuum reflection symmetric LRS solution of 
Bianchi type III of the Einstein- Vlasov equations for massless particles is rep- 
resented as a solution of ([ll]) with e = — 1 then for r — > 00 it converges to the 
the point P9 with coordinates (00, |, i) and for r — * — 00 either 

(i) it converges to Pi or 

(ii) it converges to the point P3 and it belongs to the unstable manifold of P3 
or 

(iii) it converges to P 4 

All of these cases occur, and (iii) is the generic case in the sense that it occurs 
for an open dense set of initial data. 

Proof The inequality (|lj) with e = — 1 implies that a solution of ([ll]) of Bianchi 
type III remains in a compact set and hence exists globally in r. The quantity 
B is positive in the region where B 2 < j + — 2E + ) 2 . Call this region 
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G. In the complement of the closure of G the inequality B < holds. Thus 
any stationary point with B > must lie on the boundary of G. A stationary 
point with a finite non-zero value of q must satisfy £+ = and this implies that 
£ = | , a contradiction. Thus the only stationary points occur for B = (these 
are the well-known Bianchi type I stationary points), q = or q = oo. In fact 
the only stationary points which are not of type I are those with coordinates 
(0, ^, ^) and (oo, \, \). Call these P§ and Pg respectively (see figure 3). 

The boundary of G is connected and so £ has a constant sign there. Checking 
at one point shows that this sign is positive. As a consequence, a solution can 
never leave G as r increases or enter G as r decreases. A solution which lies on 
the boundary of G at some time (with q non-zero and finite) must immediately 
enter G to the future and enter the interior of its complement to the past. 
Consider now the behaviour of a given solution as r decreases. If it stayed in G 
for ever then B would have to increase as r decreases. On the other hand, any 
a-limit point would have to be in the boundary of G due to the monotonicity 
properties of B. This is not consistent. Thus as r decreases the solution must 
reach the boundary of G and, as a consequence the interior of the complement 
of G. In the latter region B is strictly monotone and so the a-limit set must 
be contained in B = 0. Then the same analysis as in the proof of Theorem 
5.1 shows that the solution belongs to one of the cases (i)-(iii) of the theorem. 
Next consider the behaviour as t increases. As t tends to infinity the solution 
must tend to the boundary of G. If it stays in the interior of the complement 
of G then it must tend to the boundary of G as r tends to infinity and, more 
precisely, to one of the points Ps or Pg. Since E + is positive at these points, 
q — > oo and so only Pg is possible. Now suppose that the solution does meet 
the boundary of G and hence enters G itself. Then it remains in G and B is 
once again strictly monotone. As before, it can be concluded that the solution 
converges to Pg as r — > oo. 

The point Pg corresponds to a self-similar solution of the vacuum Einstein equa- 
tions much as does P-j (this time the trajectory is the horizontal straight line 
in figure 3 from Ps to Pg). This is the Bianchi III form of flat space (see p. 
193 of Jll|). Once again the nature of the initial singularity is similar to that 
in solutions of type I. On the other hand the final singularity is qualitatively 
different from any we have seen so far. In this case the solution is approximated 
at large times by a vacuum solution in the sense that the dimensionless quantity 
p/(trk) 2 tends to zero as t — > oo. A very similar analysis applies to the system 
for a radiation fluid. LRS Bianchi type III fluid solutions with equation of state 
p = behave like solutions of Bianchi type I near the initial singularity. They 
are approximated at large times by the same vacuum solution as in the case of 
kinetic theory. The approach of ||] should allow similar statements to be proved 
for other fluids with a linear equation of state, but this does not seem to have 
been worked out explicitly in the literature. 
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Figure 3: The (q, E+,2?) space and the fixed points for Bianchi type III. 
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6 Conclusions 



The above theorems show that solutions of the Einstein- Vlasov equations with 
high symmetry exhibit a wide variety of asymptotic behaviour near a singular- 
ity and in a phase of unlimited expansion. They can have a point singularity, 
barrel singularity or cigar singularity or they can show oscillatory behaviour 
near a singularity. In an expanding phase they can resemble a fluid solution 
(Bianchi type I and VIIo), a vacuum solution (Bianchi type III) or a solution 
of the Einstein equations with null dust (Bianchi type II). There are notable 
differences in comparison with a fluid model, and this includes the radiation 
fluid, which is often used as an effective model of massless particles in cosmol- 
ogy. The most striking qualitative difference is the appearance of oscillatory 
behaviour in type II solutions. It is interesting to compare this with the analy- 
sis of spacetime singularities by Belinskii, Khalatnikov and Lifschitz S. They 
do not say precisely what they assume about matter but it seems that they do 
assume, at least implicitly, that pressures cannot approach the energy density. 
This assumption is not necessarily satisfied in a kinetic description. The mean 
pressure cannot exceed one third of the energy density but if it all concentrates 
in one direction the pressure in that direction can approach the energy density. 
This leads to a source of oscillations beyond those taken into account in || . 

While the oscillatory behaviour of cosmological models near a singularity has 
often been observed numerically and explained heuristically, it has rarely been 
captured in rigorous theorems. To our knowledge the only example where this 
had been done previous to Theorem 4.2 of this paper is in a class of solutions 
of the Einstein-Maxwell equations of Bianchi type VIo analysed in . 

The results of this paper concern only massless particles. One may ask what 
would change in the results if the case of massive particles is considered. In 
one case the answer is known, namely in Bianchi type I. There the solution 
approaches a dust solution in the expanding phase. It is reasonable to expect 
that this happens more generally. As the model expands in all directions the 
pressures should become negligible with respect to the energy density, leading 
to a dust-like situation. However, the techniques necessary to prove this are not 
yet known. Near the initial singularity, the equations for massive particles look 
like those for massless particles and it may be conjectured that the behaviour 
near the singularity is similar in both cases. Unfortunately that has also not 
yet been proved. 

It is interesting to note that matter seems to have the effect of making 
the evolution of the geometry under the Einstein equations less extreme in a 
phase of unlimited expansion. In Bianchi type I the vacuum solutions (Kasner 
solutions) are such that some spatial direction is contracting or unchanging 
in the expanding time direction (the time direction in which the volume is 
increasing). This is no longer the case when perfect fluid or kinetic matter 
is added, since then the model isotropizes. In type II there is no complete 
isotropization but it is still the case that with fluid or kinetic matter all directions 
are eventually expanding, in contrast to the vacuum case. The type III case is 
borderline, since there solutions with collisionless matter are asymptotic, in the 
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sense of the variables used in this paper, to a vacuum solution in the expanding 
time direction. The vacuum solution is such that the scale factor b is time 
independent. In the other LRS Bianchi type III vacuum spacetimes this scale 
factor is asymptotically constant as t — > oo and for a radiation fluid this is 
also the case (cf. jl9| p. 203). On the other hand for dust models which also 
converge to the same vacuum model in terms of the Wainwright-Hsu variables, 
this scale factor grows without bound, although much more slowly than the 
other scale factors p. 202). It is difficult to decide what happens in 

the case of collisionless matter with massless particles, since the point Pg is a 
degenerate stationary point of the system ( |lT| ) . In the corresponding system for 
a radiation fluid the point with these coordinates is also a stationary point but 
is non-degenerate. 

For the Einstein- Vlasov equations with massless particles the LRS reflection 
symmetric solutions of Bianchi types I, II, III, VIIo and Kantowski-Sachs type 
have now been analysed as far as to give a full description of their general 
behaviour near the singularity and in a phase of unlimited expansion. There are 
still plenty of open questions related to this. What happens with LRS solutions 
of types VIII and IX? What happens if reflection symmetry is dropped? Docs 
this lead to a new kind of oscillatory behaviour? What happens if the LRS 
condition is dropped? (This is still open even in the Bianchi I case.) Can the 
Hamiltonian formulation of the equations, which played an important role at 
one point in our arguments, usefully be applied in some of these more general 
cases? Answers to these questions could help to deepen our understanding of 
the dynamics of solutions of the Einstein equations with matter in general. 
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Appendix: Some background on dynamical sys- 
tems 

First some terminology will be introduced. We use the phrase 'dynamical sys- 
tem' as a synonym for 'system of ordinary differential equations'. The difference 
between the two is then only one of point of view. A stationary point of a dy- 
namical system is a time-independent solution. An orbit of a dynamical system 
is the image of a solution. A point x* is an a-limit point of a solution x(t) 
if there is a sequence of times t n with t n — > — oo such that x(t n ) — > x*. The 
set of all a-limit points of a solution is called its a-limit set. The analogous 
notions of cj-limit point and w-limit set are obtained by replacing t by —t in 
these definitions. Basic properties are that the a-limit set is closed and that, if 
the solution remains in a compact set as t — > — oo, it is connected. If x* is a 
point of the a-limit set of an orbit then the orbit through x* lies in the a-limit 
set of the original orbit. Analogous statements hold for the ui- limit set. For 
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details and proofs see e.g. 0, chapter 7. 

If xo is a stationary point of a dynamical system we can linearize the system 
about xq. The linearized system is of the form dx/dt = Ax, for a matrix A. 
Associated to A is a direct sum decomposition E\ ® E 2 © E 3 where the vector 
spaces Ei , E 2 and E 3 are spanned by generalized eigenvectors of A correspond- 
ing to eigenvalues with positive, zero and negative real parts, respectively. These 
spaces are called the unstable, centre and stable subspaces. For each of these 
three subspaces there is a manifold which is tangent to the corresponding sub- 
space at xq and is left invariant by the dynamical system. These manifolds are 
called the unstable, centre, and stable manifolds respectively. The unstable and 
stable manifolds are unique while the centre manifold need not be. For details 
see the appendix of 0] . 

The behaviour of solutions of a dynamical system near a stationary point is 
described by the reduction theorem. 

Theorem Al (Reduction theorem) Let xq be a stationary point of a C 1 dy- 
namical system. Then the system is topologically equivalent near xq to the 
Cartesian product of a standard saddle with the restriction of the flow to any 
centre manifold. 

This theorem is proved in Topological equivalence means that there is a 
homeomorphism which takes one system to the other. A standard saddle is the 
dynamical system on R ni +™ 2 given by dy/dt = y, dz/dt = — z, where y € R m 
and z £ R™ 2 . The special case (hyperbolic case) where the centre manifold is 
trivial is the Hartman-Grobman theorem [Q . 

The next result is intuitively rather obvious, but since we do not know a 
published proof we will provide one here. 

Lemma Al Let p be a hyperbolic stationary point of a dynamical system which 
belongs to the a-limit set of a given orbit. Then either each neighbourhood of p 
contains a segment of the orbit which is contained in the unstable manifold of p, 
or the a-limit set contains a point of the stable manifold of p other than p itself. 
The analogous statement with the roles of the stable and unstable manifolds 
interchanged also holds. 

Proof By the reduction theorem we can assume that in a neighbourhood of p 
the system takes the form: 

dx/dt = x, dy/dt = -y (29) 

with solution 

x = Ae\ y = Be- f (30) 

The unstable and stable manifolds are given by y = and x — respectively. 
Suppose that there is a neighbourhood of p where there is no segment of the 
orbit contained in the unstable manifold. Then there exists a sequence of points 
p n on the orbit with non-vanishing y coordinate which converges to p. If we 
denote the coordinates of corresponding segments of the solution by (x n ,y n ), 
then y n = B n e~ l for some B n ^ 0. Consider now a coordinate closed ball 
contained in a neighbourhood of p where the reduction can be carried out. As 
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t decreases each of the solutions (x n ,y n ) must leave this ball and so must, in 
particular contain a point of the boundary sphere. Call the resulting sequence 
of points of the sphere q n . By compactness q n has a subsequence converging to 
a point q. The point q belongs to the a-limit set. Now x n = A n e* for a sequence 
with A n — > 0. Hence the x coordinate of q is zero and q belongs to the stable 
manifold of p. The proof in the case that the roles of the stable and unstable 
manifolds are interchanged is very similar, using the points where the solution 
exits the ball in the positive time direction. 

The following variant of Lemma Al allows a centre manifold of a certain type. 
Lemma A2 Let p be a stationary point of a dynamical system which belongs 
to the a-limit set of a given orbit. Suppose that the centre manifold is one- 
dimensional and that there is a punctured neighbourhood of p in the centre 
manifold which contains no stationary points and such that the solutions on the 
centre manifold approach p as t — > oo on one side of p and as t — > — oo on the 
other side. Suppose further that the stable manifold is trivial. The boundary 
between points on orbits which converge to p while staying in a small neigh- 
bourhood of p as t — > — oo and points on orbits which do not is the unstable 
manifold. The analogue of Lemma Al holds, where the stable manifold is re- 
placed by the half of the centre manifold on one side of the unstable manifold. 
This half of the centre manifold is unique. 

Proof By the reduction theorem we can assume that in a neighbourhood of p 
the system takes the form: 

dx/dt = F(x), dy/dt = y (31) 

for some function F which vanishes together with its derivative at the origin, and 
is positive otherwise. The boundary hypersurface is given by x — 0. The half of 
the centre manifold referred to in the statement of the theorem corresponds to 
x < and y — 0. In the half-plane x < the system is topologically equivalent 
to a hyperbolic saddle and so it is possible to obtain the conclusion as in the 
proof of Lemma Al. 

Next we state the Poincare-Bendixson theorem. The form of this theorem 
which we will use is the following (cf. [7], p. 151): 

Theorem A2 (Poincare-Bendixson) Let U be an open subset of R 2 and consider 
a dynamical system on U without stationary points. Let x(t) be a solution which 
exists globally and remains in a compact subset of U as t — > — oo. Then the 
a-limit set of the given solution is a periodic orbit. 

The analogous statement holds for the the w-limit set. A periodic orbit is, of 
course, just the image of a periodic solution. 
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